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ABSTRACT 

Owing to their more extensive sky coverage and tighter control on systematic errors, future deep weak 
lensing surveys should provide a better statistical picture of the dark matter clustering beyond the level 
of the power spectrum. In this context, the study of non-Gaussianity induced by gravity can help tighten 
constraints on the background cosmology by breaking parameter degeneracies, as well as throwing light 
on the nature of dark matter, dark energy or alternative gravity theories. Analysis of the shear or flexion 
properties of such maps is more complicated than the simpler case of the convergence due to the spinorial 
nature of the fields involved. Here we develop analytical tools for the study of higher-order statistics such 
as the bispectrum (or trispectrum) directly using such maps at different source redshift. The statistics we 
introduce can be constructed from cumulants of the shear or flexions, involving the cross-correlation of 
squared and cubic maps at different redshifts. Typically, the low signal-to-noise ratio prevents recovery 
of the bispectrum or trispectrum mode by mode. We define power spectra associated with each multi- 
spectra which compresses some of the available information of higher order multispectra. We show how 
these can be recovered from a noisy observational data even in the presence of arbitrary mask, which 
introduces mixing between Electric (E-type) and Magnetic (B-type) polarization, in an unbiased way. 
We also introduce higher order cross-correlators which can cross-correlate lensing shear with different 
tracers of large scale structures. 
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> INTRODUCTION 



r eak lensing surveys play an important role as cosmological probes complementary to Cosmic Microwave Background (CMB) surveys and large-scale 
laxy sur veys. In principle, th ey can probe evolution of the dark matter power spectrum at a moderate redshifts i n an unbiased way; for a recent 
•review, seelMunshi et al.l d2008h. However, the first weak lensing measurem ents were published within the last decade (Bacon, Refregier & Ellis 2000; 
Wittman et alll2000l ; iKaiser. Wilson & Luppino 2000; Waer beke et al|[2000h so this is a young field. There has been a tremendous progress in the last 
decade since the first measurements, on various fronts, including analytical modeling, technical specification and the control of systematics. Weak lensing 
probes cosmological perturbations on smaller angular scales where the perturbations are both nonlinear and non-Gaussian. The two-point correlation 
function (or, equivalently, the power spectrum) of weak lensing is weakly sensitive to the cosmological constant ..a- It also depends only on a degenerate 
combination of amplitude of matter power spectrum erg and the matter density par ameter Qm- This degeneracy can be broken by use of the three-point 
correlation function which is diagnostic of the non-Gaussianity induced by gravity jVillumseiJ 19961 ; H am &Seliakll997h . 

The higher order statistics needed to analyze weak lensing in detail are however prone to the effects of noise due to the intrinsic ellipticity distribution, 
finite number of galaxies (shot noise) and sample (or cosmic) variance owing to the partial sky coverage. Ongoing and planned weak lensing surveys, 
such as the CFHT leg acy surve£l Pan-STARRsH the Dark Energy Survey, and further in the future, the Large Synoptic Survey Telescope^, IDEM and 
Euclid will provide a wealth of information in terms of mapping the distribution of mass and energy in the universe. A large fractional sky coverage as 
well as accurate p hotometric redsh ift determination will help to probe higher order correlation functions at a higher signal-to-noise(S/N) then is possible 
at present (see e.g. lPen et alj d2003tv ). 

Almost all previous studies concerning weak lensing have used a real space description and employed correlation function of various orders to probe 
underlying mass clustering. Future surveys will cover a good fraction of the sky, if not the entire sky. This motivates us to use a harmonic description and 
employ the mutispectra - which are simply the Fourier representations of the higher order correlation functions. Weak lensing observations require masks 



1 http://www.cfht.hawai.edu/Sciences/CFHLS/ 

2 http://pan-staiTS.ifa.hawai.edu/ 

li http://www.lsst.org/llst_home.shtml 
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with complicated topology to deal with presence of foreground objects. Our formalism in harmonic space can deal with arbitrary mask. The formalism 
described here is completely general and can deal with fields with arbitrary spins including shear and flexions. 

In the absence of information concerning the redshifts of the source galaxies, traditional weak lensing studies have tended to operate in p rojection 
or in 2D slices. Indeed, before the advent of 3D weak lensing, almost all studies were carried out in such a w ay (Jain, Seljak & White 2000). Due to 
lack of all-sky coverage most studies also used a "flat sky" approach dMunshi & Jainl (BoOll); iMunshil ( feOOOh: Munshi & Jain ( 200 CQ) which has le d 
to study of lower order non-Gauss i anity using convergence k as well a s shear "/ IMunshi & Jainl d20o1E : IValageasT (2000): iMunshi & Valageasl (2005): 
IValageas. Barber. & MunshU d2004h : lvalageas. Munshi. & Barbell d2005l) . As the next sta ge of developmen t tomo graphic approac hes went beyond pro- 
jected survey by binning galaxies in redshift slices which can further tighten the constraints lTaka d a & W hite (2003); Tak ada & Jainl J20041) . More recently, 
the use of pho t ometr ic redshifts to study w e ak lensing in three dimensions was i ntroduced b y Heavens (2003). It was later developed by many authors 
iHeavens et al. (2000); Heavens et all d2006l): iHeayens. Kitching & Verdd d2007l): ICastro et all d2005l) . and was shown to be a vital tool in constraining 
dark energy equation of state heavens et alll2006l) . neutrino mass dKitching et al.ll2008h and many other possibilities. In our present study we propose 
estimators for non-Gaussianity using projected data. 

In most cosmological studies, the power spectrum remains the most commonly used statistical probe and its evolution and characterization remain 
the best studied. As has already b een pointed out, however, higher order statisti c s can help to break degeneracies e . g. in Qm and ag, but a r e more difficult 
to prob e observationally (see e.g. lBernardeau. van Waerbeke & Mellied ( fl997h : |jain & Seliakl dl997h : lHull d 19991) : Schne ider et a3 d2002l) ; lT7kada& Jainl 
(2003)) given the low signal to noise associated with them. As a result most com monly used high er order probes compress the available information 
to a simple number (e.g. "skewness" or "kurtosis") which are one-point statistics dPen et alj|2003l). The higher order correlation functions have been 
detected observationally dBernardeau. van Waerbeke & Mellierll997l:lBernardeau. Mellier & Van Waerbeke 2002) and future larger samples are expected 
to improve the statistical situation still further. A recent study bv lMunshi. Heavens & Colesl 1 201ol) has underlined the usefulness of working with two- 
point correlations (or their Fourier analog the power spectrum) at each order. At the level of the bispectrum, a single power spectrum tends to compress 
some of t he available information; there is m ore than one degenerate power spectrum associated with the higher order multispectra. We extend the recent 
study by IMunshi. Heavens & Colesl d2010h and provide systematic results at each order for spin— 2 fields, thus generalizing our previous results for 
spin-0 fields s uch as the convergenc e. Understanding higher order statistics also provi de an ind i cation of t he scatter associate d with estimation of lower 
order statistics iTakada & Jainl d2009l) . Their resu lts extend the fo rmalism developed in Heavens (2003f) and lCastro et all d2005l) and we follow the all-sky 
expansion introduced in weak lensing studies bv lStebbinsl ( [T993) . 

Statistics of the shear 7 or convergence k seen in a lensing map are sensitive to the statistics of underlying density contrast S — 5pb / pt of background 
density pb. An accurate understanding of the gravitational clustering is therefore essential for modeling the weak lensing statistics. At present we lack 
a detailed analytical model of the non-linear gravitational clustering. In the absence of a clear a n alytical pic t ure, b e yond the perturbative reg i me th e 
modeling of the higher order s tatistics are done u si ng a hierarchical ansatz (see e.g. iFrvl dl984l); ISchaeffed dl984l): iBernardeau & Schaefferl dl992h ; 



Szapudi & Szalavldl993l . [l997l) : lMunshi et all dl 9991) : lMunshi, Coles & Melottldl999jbl) ; IMunshi, Melott & Colesl dl999l) ; lMunshi & Colesl d2000L[200l 



2003)). The hierarchical ansatz assumes a factorisable model for the higher order correlation functions. Different hierarchical ansatze differs the way they 



assign amplitudes to the various tree amplitudes. We will employ a very generic form of the hierarchical ansatz which has already been tested in modeling 
of projected statistics. The method is flexible enough to allow for a more specific prescription. Other approaches to model nonlinear gravity include the 
halo models dCoorav & Seth 200 j) which has also been employed in weak lensing studies. 

In a recent paper, IMunshi. Heavens & Colesl d2010h advocated the use of higher order cumulants and their correlators of convergence k along with 
their associated power spectra for studying dark matter clustering beyond power spectra in 3D. The aim of the present paper is to extend that study to 
fields with non-zero spin (e.g. shear). Shear can directly be used from observed data to study the dark matter clustering beyond the usual approximation of 
Gaussianity. Although in this study we primarily focus on shear, the analytical results are directly applicable for hig her order spin fields sometime used in 
weak lensing known as flexions. Previous studies have also focused on the cumulant correlators of the shear 7 and k ( Bernard eau. van Waerbeke & Mellierl 
ll997l ;rBernardea u. Mellier & Van Wae rbeke 2002; Bernardeau, Mellier & van Waerbeke 2003). While cumulant correlators for convergence fields can be 
probed in a relatively straight forward manner, even to arbitrary order dMunshil2000l) . the analogous studies are more involved for she ar fields. There have 
been some studies involving smoothed M ap statistics to study cumulant correlators from shear map by using compensated filters ( Munshi & Valageas 
2005), but in this present paper we develop a more systematic approach to probe higher order statistics directly through spinorial objects. The focus of the 
present study is to develop tools in the Fourier domain and to use the power spectra associated with the multispectra as opposed to cumulant correlators. 
Given the non-trivial topology of weak lensing surveys correlation functions indeed are a good option to work with, but, as we will show, the power 
spectra we develop can also be estimated in the presence of arbitrary mask and noisy data. 

The paper is arranged as follows. In §2 we discuss the basic formalism and introduce some terminology and notation. In §3, we introduce the 
power spectra related to bispectra and trispectra for various combinations of the convergence, shear as well as flexion fields. In §4 we focus on statistical 
description of underlying matter distribution and relate them to power spectra defined in previous sections for weak lensing observables. Finally §5 is 
reserved for discussions and future prospects. 



2 CONVERGENCE, SHEAR AND FLEXIONS 

Among the observables including in weak lensing studies are the convergence shear ~y(£l) and flexions J-(Cl), Q(£l). The harmonic decomposition 

needed in a l l-sky calculations involves spin-weight spherical harmonics for such objects. We introduce the basic notations in this section following 
ICastro et H d2005l) . We will ignore the radial dependence for this work as we mainly concentrate on projected surveys. The results presented here 
will extend the results in dMunshi. Heavens & Co les 2010) wherein only scalar fields were considered. Though the analysis here follows full spherical 
decomposition it is indeed possible to choose a rectangular coordinate system which can simplify the calculation for surveys observing a small patch of 
the sky. 
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Power Spectrum 




Figure 1. The plots correspond to the power spectrum of convergence C ; KK . The source redshift is z\ = 22 = 1. See text for more details. 



2.1 Spin-weight Spherical Harmonics 

We start by introducing the spin-we i ght s p herical harmonics aYi m (Cl) th at can be expressed in terms of the Wigner-D function (see e.g. 
IVarshalovich. Moskalev & KhersonskiH Jl988h ; IPenrose & Rindleil fl984,1986) for a detailed discussion). They form a complete and orthonormal set 
of function on the surface of the celestial sphere. Any spin-weighted function s /(") can be decomposed into harmonics s fim using spin-weight spherical 
harmonics: 



s y im (n) = J^^DL. im (d,^,oy, j(Ci) = ^ s Y lm (£i) a f lm . (i) 

lm 

The spin- weighted sp herical harmonics s Yi m (Cl) (define d only for \s\ > I )can be seen as a generalization of the scalar, vector and tensor, spherical 
harmonics ( Varshalovich, Moskalev & Khersonskii 1988) satisfying the following orthogonality and completeness relations: 

^2 s Y lm (h) a ,Y lm (n') = s D (Ci - n'); J dh .Y, m (ft).,y, w (fi) = sf s ,s§6^ mt . (2) 

lm 

Here 8 K is the Kronecker delta function and 8d is a Dirac delta function. The spin raising operator "edth", 9 and its complex conjugate, the spin lowering 
operator 8 are defined through their action on spin spherical harmonics which can be used to define spin-weighted harmonics 3 Yi m (fi) from ordinary 
spherical harmonics Yi m (Cl). The spin spherical harmonics s Yi m (£i) are the eigen functions of the operator 9 8 . Note that 9 can be thought also as an 
effective covariant derivative on the surface of the sphere which generates quantities related to the various spins: 

QsY lm = [{I - s)(l + s + 1)] 1/2 s+1 Y lm ; B s Y im = [(l + s){l - s + 1)] 1/2 s -xYi m . (3) 

This also implies that 9 9 a Yi m = [(7 — s)(l + s + l)] s Yi m and 9 9 s Yi m = [(/ + s) (i — s + l)] s Y! m - The conversion between convergence and shear 
harmonics used the following expression which we will also be using later: 



g gy im (n) = 




The operators 9 and 9 are convenient for expressing the spin weight objects which are coordinate frame dependent in terms of ones which are frame 
independent or scalars. All invariant differential operators on the sphere may be expressed in terms of 9 and 9 . 

The following overlap integrals will be useful in our future derivations which will be used to derive the expressions for bispectrum and trispectrum 
related power spectrum. 

( s Y lm {Cl) s ,Y Vml {tl) a , l Y l u m , l {(l)dh = S wvl ( 1 l ', l "„)( 1 l ' , l "„); s„n>'=\/¥^ ){2l ' + mi '' + 1) ~. (5) 
J \ m m m J \ ~s ~s — s J V 4-7T 

The above result can be cast in following form which is useful for expressing integrals of more than three spin spherical harmonics in terms of Wigner 3j 
symbols when dealing with trispectrum or multispectra of even higher order. 
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s Y lm ^)MCl) = £ sYl M S Llll „ f L M l m , ) ( L s J' s , l " ) . (6) 

LSM ^ ' ^ ' 

These results are generalizations of similar results often used in the context of scalar (spin-0) harmonics which can be recovered by setting s — s' = 
s" — 0. In a small patch of sky, a rectangular coordinate system can be introduced which can simplify some of the calculations and might be useful for 
surveys with small sky coverage. A detailed analysis will be presented elsewhere (Munshi et al. (2010) in preparation). 



2.2 Harmonic decomposition of Convergence and Shear 

We start by defining the complex shear 7 in terms of the individual shear components 71 and 72 at a given angular position Q = (8, <f>) an d at a radial 
line-o f- sight distance | r | = r (see Mu nshi et alj J2008h for a more detailed discussion and definitions). We will closely follow the notations of ( Cas tro et"a3 
l2005h where possible. We will not Fourier decompose in the radial direction. The harmonic decomposition is performed only the surface of the celestial 
sphere 



7 ±(A,r)=7 1 (n,r)±iTia(fl,r). (7) 

It is customary to define a complex lensing potential 4>{r) = 4>e(t) + icj>B (r) which is related to the line-of-sight integration of the peculiar gravitational 
potential <j> and to define the shear in terms of spin-derivatives 9 and its conjugate 9 of <j>. The action of is to r a ise the spin whereas 9 reduces t he spin 
of the object it is acting on. This operator was first introduced by Varshalovich, Moskalev & Khersonskiil il988h ; |Penrose & Rindled {1984,1986) on the 
surface of the sphere in order to define the now widely used spin-weight s spherical harmonics which live on the surface of a sphere 

7 (r) = iag[^(r) + ^ B (r)]; 7 *M = ~88[<^(r) - #s(r)]. (8) 

Explicit expressions for these spin-lowering and raising differential operators are given in lCastro et all {2005). While acting on a complex scalar potential 
<f>, the fields 7(r) generates the 7(r) and its conjugate which are of spin +2 and —2 respectively. In later sections we will use the generalized symbol S T 
for general spin fields which will include products of shear fields as well as higher derivative spin fields such as flexions. In our current notation 2V — 7 

and _2F = 7*: 

2 r(r) = 7 (r) = iga^(r); _ 2 r(r) = 7 *(r) = ^880* (r). (9) 

The individual shear components 71 (r) and 72(1") can be expressed in terms of a complex lensing potential <f>(r) = 4>e(t) + i</>s(r). The magnetic part 
of the potential (f>B (r) will take contribution mainly from systematics and electric part corresponds to pure lensing contribution 

7 i(r) = i(gg + 86Mr); <ys« = -J(3 3 - 8 S)#r); «(r) = ~( 9 5 + 9 5 )0(r). (10) 
The harmonics of 7 and 7* can be expressed in terms of the harmonics coefficients of $_b and $_b denoted by Ei m and Bi m respectively. 

2Tl m — —[Elm + iBi m ]; -2T; m = — [Ei m — iBlm]. (11) 

The spinorial fileds can likewise be expanded in an appropriate basis which uses the spin-spherical harmonics s Yj m (fl). For the case of 2 r; m we have 
the following expression: 



±2r(r) = X) E \l 773§T ±* Y l™(r, ft). (12) 

1=0 m=-l V ^ 

We will be only dealing with 2D part of the spherical expansion as we will be dealing with projected surveys in this work. The harmonics Ei m and Bi m 
are harmonic components of Electric (E) and Magnetic (B) fields respectively. In the absence of B-mode the harmonics of the shear components are 
directly related to the harmonic component of the Electric field Ei m . The harmonic transforms of the shear components and convergence are linked to 
the lensing potential (f> as follows: 



These expression are useful in relating the bi- or trispectrum of various quantities in terms of the convergence ones. 
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2.3 Harmonic Decomposition of Flexions 

Highe r spin objects known as flexions are sometimes al s o used to study weak lensi n g; they are related to de rivatives of the shear dGoldberg & NataraianI 
d2002h : lGoldberg & Beacon! d2005l) ; lBacon et all d2006h : lBacon & Goldberg! d2005l) : ISchneider & End2008h ). There are two flexions which are typically 
used; T also known as the first flexion (spin -1) and Q which is also known as the second flexion (spin -3). The combinations of these two flexions can 
specify the weak "arciness" of the lensed image an d hence can quantify the lensing distortion beyond that encoded in sh ear. Their relationship with the 
shapelet formalism have been discussed at length (Refregier 2003; Bernstein & Jarvis 2002; Refregier & Beacon 2003). Both flexions have been used 
extensively in the literature for individual halo profiles and also for the study of substructures (Bacon et aL 2006). Our aim here, as it is in the case of shear, 
is to focus mainly on higher order statistics of these objects for generic underlying cosmological clustering. We will do so by expressing the bispectrum 
of the shapelet s in terms of that o f the convergence k, as we did for shear 7. The first and the second flexions T and Q can be derived from the lensing 
potential (j>(r) dCastro e t al 2005) by the use of the lowering and raising operators: 

J"(r) = ~ (9 9 9 + 9 5 9 + 9 9 B) 4>(r); G(r) = i 9 9 9 </>(r). (14) 

Flexions have been used primarily to measure the galaxy-galaxy lensing to probe the galaxy halo density profiles. Their cosmological use will depend on 
an accurate understanding of gravitational clustering at small angular scales. In Fourier space we will denote the harmonics of T and Q by Tim and Qi m 
and we can use the Eq.10 to express them in terms of the 4>i m . 



Ti m = ±l 1/2 (l + l) 1/2 (3l 2 +3l-2)<j }lm ; G lm = ij|±|ll0 im (15) 

Though we will primarily be focusing on higher order statistics of shear the results are equally applicable for flexions (which generalize the concept of 
shear). In addition to the shear statistics we will also consider a scalar field in our analysis We will leave this, which can represent a suitable large scale 
tracer, arbitrary. The resulting analysis depends on the multispectra involving both the weak lensing shear field and that of the tracer fields. The flexions 
put more weight on smaller scales, where the hierarchical ansatz is known to be more accurate, so the formalism we develop here will also be suitable for 
them. It is also possible to use our general results to construct estimators for mixed bispectra involving shear and flexions. 

However we should also point out that the modeling of intrinsic flexions of source galaxies, which is the main contribution to the noise, is difficult. 
This will depend heavily on modeling of galaxy shapes beyond the simplest description of their intrinsic shear. This uncertainty is expected to increase 
with survey depth. A detailed modeling and optimization of survey strategy will be presented elsewhere. 



3 HIGHER ORDER STATISTICS OF WEAK LENSING SHEAR AND FLEXION 

The Higher order statistics encodes the departure from Gaussianity. The non-Gaussianity probed by any cosmological observations can either be a primary 
one due to initial conditions or a secondary one which in case of weak lensing is mainly induced by subsequent gravitational instability responsible for 
structure formation. At lower redshift it is expected that the secondary non-Gaussianity will dominate the primary one. The estimation of non-Gaussianity 
using higher order correlation functions are typically dominated by noise. The higher order correlation functions or their Fourier transforms (multispectra) 
contains a wealth of information. Some of these information however can be degenerate due to large number of possible configurations in real or Fourier 
domain. Typically one-point statistics e.g. skewness or kurtosis that are used compress all possible information to a single number thereby increasing 
the Signal-to-Noise b ut reducing th e information content. As a next step, collapsed two-point objects (or cumulant correlators) which are higher order 
objects are also used (Munshi 2000). They are easier to recover from a noisy data than the complete correlation functions. The Fourier transforms of these 
two-point cumulant correlators are the power spectrum associated with a specific choice of the multispectra. There can be more than one power spectra 
for any specific choice of multispectra. In this section we derive power spectra associated with bispectra and trispectra for shear as well as flexions. 

3.1 Power spectrum associated with the Bispectrum 

The bispectrum is the lowest in the hierarchy of multispectra and have the highest Signal-to-Noise. There has been lot of interest in modeling and 
detection of weak lensing bispectrum mainly using convergence maps as well as shear. As remarked earlier, shear correlation functions are more 
complicated and have rich patterns given t heir spinorial nature. They have been stu died in real space for specific triplet confi gurations analytically 
dBernardeau. Mellier & Van Waerbeka 2002) and have b een measured in re al surveys dBerna rdeau. M ellier & van Wae rbeke 2003) f° r fi rst detection of 
non-Gaussianity signal from weak lensing observations. iBernardeaul d2005h has also provided analytical results for more general configurations but in 
an idealistic one-halo configuration. In this paper we focus on the power spectra associated with bispectra to probe higher order correlation hierarchy in 
studying gravity induced non-Gaussianity. Our results are completely general and given any models of underlying bispectrum can predict observed power 
spectrum related to bispectrum in the presence of arbitrary observational mask. We will assume noise to be Gaussin through out and will not contribute 
to the bispectrum. 

3.1.1 The Two-to-One Power spectrum: C I ' 

We will start by considering two fields on a surface of a sphere s r(fl) and s iT'(Cl), which are respectively of spin s and s'. These objects can be the 
shear, j±, flexions, T, Q or a tracer field $ '. The resulting product fields such as e.g. 7+7- can be of spin zero, on the other hand 7+ and 7I are of 
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Bispectrum 
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Figure 2. The plots correspond to the bispectrum related power spectrum (7; ■ Various curves correspond to the spectra Cf C^ K ' 1+ , C^ 1+ ,7+ and c 7+7 +' 7 + as 
depicted. The source redshift is z\ = Z2 = 1. See text for details. 



spin —4 and +4 respectively. The results can also involve fields such as which is a spin —2 field. A spin s field can be decomposed using spin- 
harmonics s Yi m {Q,). These spin harmonics are generalization of ordinary spherical harmonics Y; m (f2) which are used to decompose the spin-0 (or scalar) 
functions e.g. defined over a surface of a sphere (sky). Throughout we will be using lower case symbols s, s' to denote the spins associated with the 

corresponding fields denoted in italics. The fields that are constructed from various powers of 7+ and 7_ can be expanded in terms of the corresponding 
spin-harmonics basis s Yj m (f2). The product field such as [ s r(r) s ;r'(r)] which is of spin s + s' can therefore be expanded in terms of the harmonics 
s+s'Yi m (£l). The following relationship therefore expresses the harmonics of the product field in terms of the harmonics of individual fieldsQ. 

[,r(r).,r'(r)] hn = j dn s r(r) s /r'(r)[ s+s ,y4(Q)] = ^ s r iimi (r) s ,r; 2m2 (r) J [.yi imi (A)][*yj araa (fi)]U^(6)]dft; 

V- r / s w / w ( h h I \( h h I \ T / (2Z! + 1)(2Z 2 + 1)(2Z + 1) 

= ^ s r iiroi (r) s ,r i2m2 (r)/ !li2! [ g/ _ {s + sl) )l mi m2 ); /w = V ^ (16) 

The expression derived above is valid for all-sky coverage. Mixing of Electric-i5 and Magnetic-i3 modes results from partial sky coverage. The results 
will be generalized later to take into account arbitrary partial sky coverage and we will derive exact form for the mixing matrix. We have assumed a 
Gaussian noise resulting from intrinsic ellipticity distribution of galaxies, however this will not contribute to this non-Gaussianity statistic. We define the 
associated power spectra to write: 



(n,r 2 ) ^EU^VMv-Cir) (17) 



c rr',r" 



(s + s') J ' 1/(1121 V (2Z + 1 



The power spectrum C ; rr ,r (n, T2) essentially cross-correlates the product maps at a given radial distance r\ to another map at another redshift to 
probe the third order statistics bispectrum B. Similar statistics in the coordinate space has been reported before. dBernardeau. van Waerbeke & Mellieil 
120031) studied {j 2 (£l)j(£l')) which directly deals with shear maps as opposed to convergence maps. This statistics along with a similar but simpler version 
which uses was studied. They employed perturbation theory to model underlying mass distribution and used a flat sky approximation to 

simplify their calculations. A complementary statistics ((7(^1) -7(r22))7) was also considered which relies on more detailed modeling of the bispectrum. 
These statistics were used by dBernardeau. Mellier & Van Waerbekell2002l) later to detect non-Gaussianity from VIRMOS-DESCART Lensing Survey. 

Our results presented here deal with power spectrum associated with the higher order multispectra and are derived using generic all-sky treatment 
and can handle decomposition into Electric and Magnetic components in a much more straightforward manner. Our results are complementary to their 
analysis. We will also generalize these results to higher order and show how to take into account the mask in a generic way. The results presented here 
are not only applicable to shear or convergence but are also applicable for higher order spinorials such as Flexions. 

We have introduced the bispectrum B^j 2 f (ri,r2,rs) in Eq . (1 1 8ft which can now be related to the harmonics of the relevant fields by the following 
equation: 



4 It is possible to prefix the spin s objects as s r(r), however to simplify notations we will associate spin s with T(r), s' with r'(r) and so on, through out the paper, and 
the spin s should be obvious from the context. 
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Bfr'X (r u r 2 ,r 3 ) = £ { s T llTni (n).^™, Ms-H^ fo)>c ( ^ ^ M. (19) 

The matrices denote Wigner 3j symbols jEdmondsl l968) which are only non zero when the quan tum numbers Z» and rrij satisfy certain conditions. For 
x = y — this generalizes the results obtained in dCooray|2005F) valid for the case of scalars (see also lChen & SzapudU J2007h for related discussions). The 
result derived above is valid for a general E and B type polarization field. Contribution from (B-type) magnetic polarization is believed to be considerably 
smaller compared to the (i?-type) electric polarization. The results derived above will simplify considerably if we ignore the B type polarization field 
in our analysis. In general the power spectra described above will not be real though real and imaginary parts can be separated by considering different 
components of the bispectrum. However if we assume that the magnetic part of the polarization is zero the following equalities will hold: 



B JXh = F h B i*2h'-> B tit2h = F ii F h B hlfl 3 ; Bl^htf* = F h F l» F h B hhh- ( 2 °) 

The result presented above is valid for all-sky coverage. Clearly from practical considerations we need to add a foreground mask. If we consider the 
masked harmonics associated with the mask, for a given w(Cl) arbitrary mask , then expanding the product field, in the presence of mask, we can write: 

[ s T(r) s ,r'(r)w(Cl)] lm = sr hmi (r)s>r' hm2 (r)w lama J [ s Y hmi (Cl)][ s ,Y l2m2 (Q)][Y lama ][ s+s ,Y* m ]dCl 

lirrii ;Z a m n 

E E i- 1 ) 1 ' ^ r il^t( r )s'^l2m2ir)w lam Jl lhl f ^ ^ _// +fl /S ) ( ^ J m l ) /[-K^il'm'RmJs+s'C] 

J. ■ m - I _ m . _ >■ '' >• / / \ /J 



liTrii l a m 

= 12(- 1 ) l+l E ar h™i( r )s' r hm2(.' r )' w l*™* I hhl' I l'la 



h h V \f h h I' 

S s' — (s + s') J \ 771i 7T72 —17l' 



x ( V L 1 \( V L 1 ) (21) 

\(s + s') -(s + s')J\m' m a -m J ' ( ' 

In simplifying the relations derived in this section we have used the relationship Eq.lO and Eq.©. We have also used the fact that a Yi m = 
( — l) m+s Yi,- m , where * denotes the complex conjugate. Similarly we can express the pseudo-harmonics of the field Z(Q) observed with the same mask 
in terms of its all-sky harmonics. 



[ s „r"w] lm (r) = I dti[ s „r"(Ci)w(£i)][ s „Y l i m (ti)} = J2 s" r " 3 m 3 (r)w hmb J [ s „y !3m3 (n)][y !i , ro6 (n)][ s „Yi^(fi)]dQ 

liTTli 

= E E ^hm2(r)w hmb I W ( J, h Q _[„ ) ( ^ % b _' m ) . (22) 

l 3 m 3 l b m b 

The harmonics of the composite field [r(n)r'(n)] when constructed in a partial sky are also a function of the harmonics of the mask used wi m - 
The simplest example of a mask would be w = 1 within the observed part of the sky and w = outside. More complicated mask The harmonics Wi m 
are constructed out of spherical harmonics transforms. We also need to apply the mask to the third field which we will be using in our construction of 
bispectrum related power spectrum or skew-spectrum associated with these three fields. The masks in each case are the same, however the results could 
be very easily generalized for two different masks. The pseudo power spectrum C\ is constructed from the masked harmonics of relevant fields which is 
a cross-correlation power spectrum: 



1 ' 

cf r - r (n,r 2 ) = [.r(n)^r , (n)iu(n)] !m (n)[.«r"(n)t l ;(n)]L(ra). (23) 

m — — I 

The pseudo power spectrum Ci is a linear combination of its all-sky counterpart C\. In this sense masking introduces a coupling of modes of various order 
which is absent in case of all-sky coverage. The matrix M u i encodes the information regarding the mode-mode coupling and depends on power spectrum 
of the mask w\. For surveys with fractional sky coverage the matrix is not invertible which signifies the loss of information due to masking. Therefore a 
binning of the pseudo skew- or kurt-spectrum may be necessary before it can be inverted which leads to recovery of (unbiased) binned spectrum. 



C ; rr '' r "(n,r 2 ) =Y J Ml°:> e "cf> v "{rx,r*) (24) 
V 

The all-sky power spectrum Cf r ,r can be recovered by inverting the equation which is related to the accompanying bispectrum by the Eq.d 1 8b which 
we discussed before. The mode-mode coupling matrix depends not only on the power spectrum of the mask but also of the spin associated with 
various fields which are being probed in construction of skew-spectrum or the power-spectrum related to the bispectrum. 



8 Munshi et al. 



m^" = ^E(2/' + i)(x + d( s ; s , l ; o -^)w a (25) 

The expression reduces to that of temperature bispectrum if we set all the spins to be zero x — y = z = Oin which case the coupling of various 
modes due to partial sky coverage only depends on the power spectrum of the mask. In case of temperature (spin-0) analysis we have seen that the 
mode-mode coupling matrix do not depend on the order of the statistics. The skew-spectrum or the power spectrum related to bispectrum, as well as its 
higher order counterparts such as power spectrum related to tri-spectrum can all have the same mode-mode coupling matrix in the presence of partial 
sky coverage. However this is not the case for PS related to the polarization multispectra. It depends on the spin of various fields used to construct the 
bispectrum. It is customary to define a single number associated with each of these bispectrum. The skewness is a weighted sum of the power spectrum 
related to the bispectrum S^ r ' r = ^2 t {2l + l)Cf r ,r . 

We list the specific cases of interest below. The relatio ns between the bispectra and the associated power spectra generalizes previously obtained 
cases where only the temperature bisp ectrum was considered dCoora\l2005l) . These results were obtained in the context of cosmic microwave background 
(CMB) polarization studies recently dMunshi. Heavens & Colesll2010f) . The radial dependence in case of weak lensing observables needs to be properly 
taken into account. 



C l ± (ri,r 2 )=}^Ji 1 i 2 iB lll2l (r 1 ,r 1 ,r2){ ) ( 26 ) 



Mr = ix>+i)(x+i)(j ^ <• <' 2 ) KJ *. (27) 

Next we can express the power spectrum Cf E,E probing the mixed bispectrum Bf x f E by the following relation: 

C? ±M ( fl ,r J )=E4 W C(n,ri,r J )(j 2 h Q 1 ) (28) 
hh V ' 

M^ = ip^l)(^ + l)(i 2 \ 4)(J 2 \ 4)KJ 2 . (29) 
Similarly for the case (7 7 ± 7 ± ,7 ± which probes the bispectrum B E f 2 f 'we can have the following expressions. 

C? ±7± ' 7± (n,r 2 ) = E^ iy C(n,n, r2 )f £ 2 l l 2 ^ (30) 



hh 



^i£(^ + i)(x + i)( ± ' 4 \ 4)(i 2 \ 4)kj 2 - (3d 



la 



Other power spectra such as, £< 7 ± 7 ±' 7 t an( j ^ 7 ± 7 ±' can a i so ^ e (j ea ][ w ^ m m a similar manner. The simp ler case of con vergence, i.e. C™' K which 
being a spin field can now be analyzed as a special case and will reduce to the results previously obtained by JCooray||200ll . The above results display 
expressions involving shear 7 and k, similar results for flexions T and Q can be derived using the same general result Eq.lll8b. We list below the results 
involving only T or Q flexions. It is of course possible to consider mixed bispectrum involving different combinations of T and Q. The bispectra B T:Fjr 
and B sss given below can be expressed in terms of convergence bispectra using the expression Eq.J15t. i.e. we have Bff? — F^F^F^Bil^ with 
Ff defined similar to the case involving shear. 

Cr' T (ri,r2) = ^2J h i 2 iBf 1 fJ(r 1 ,r 1 ,r 2 )( l l \ [A (32) 
hh V ' 

r = ^(2l' + m i a + l)f[ \ M(J l ; l '\\ WlA \ (33) 



M x 

The power spectrum related to the second flexion Q can likewise be written in terms of its bispectrum B sss 



Cf^ S {r 1 ,T 2 ) = Y. J ^ B hZ^r 1 ,r 2 )( ^ h 3 J] (34) 



In 



4tt^^ " v a 1 _g _ 3 6 1 1 _ 3 3 



As is clear that the coupling matrix M for different cases, depends not only on the power spectrum of the mask but also on the spins associated 
with the fields that are being probed. The inversion of coupling matrix M u i can require binning. The binned coupling matrix M bb r in case of small sky 
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coverage will lead to recovery of binned power spectra C'i b . The mask is important even for small sky coverage as the observed region can have non-trivial 
topology because of presence of very bright stars or other similar foreground objects. 

The different power spectra such as C* K ' K ' and C* 77 ' 7 though intrinsically probe the same bispectrum B EEE weights individual modes in a different 
manner and can be used to probe effect of systematics in handling real data. Each of these power spectra can be used to define corresponding one-point 
skewness. We plot the results of our numerical calculations in Fig-|2]for various bispectrum related power spectrum. The power spectrum for convergence 
for the same models are displayed in Fig-Q] For these calculations we put all the sources at a same redshift z s = 1. The results can take into account for a 
source redshift distribution with a given spread and median redshift of sources. 



3.2 Power spectrum associated with the Trispecrum 

Where as the bispectrum represents the lowest order deviation from the Gaussianity, the next higher order representation is Trispectrum. Reason to go 
beyond lowest order non-Gaussianity is multifold. While lowest order non-Gaussianity indeed probed directly by bispectrum with a higher signal-to-noise, 
the study of trispectrum can go beyond lowest order in probing the gravity induced non-Gaussianity. 

(2 2) 

In a different context, relating to CMB studies involving secondaries it is known that trispectrum related power spectra Cy ' can separa te out lensing 
induced non-linearity from other contributions from secondaries without involving cross-correlations with any external data sets see e.g. l lMunshi et all 
|2010|) for extensive discussions. It can also put independent constraints on non-Gaussianity parameters g^L and /nl- The power sp ectrum Cp' 1 ^ 
associated with the Trispectr um has also been used in studies of 21cm surveys which maps neutral hydrogen distribution at higher redshift dCooravll2006l : 
ICoorav. Li & Melchiorrill2003l) . 

We explore here the possibility of using these power spectra for studies in the context of weak lensing. To model weak lensing tri-spectra one needs 
to model the underlying mass trispectra. This we do here using a well motivated model, namely hierarchical ansatz which is known to be valid at smaller 
angular scales. However we want to stress that modeling of underlying trispectra can always follow more complicated prescriptions from (extensions of) 
perturbation theory or various improvements halo models which are also employed frequently in such studies. 



3.2.1 The Two-to-Two Power spectrum: G\ 

In constructing the trispectrum related power spectrum we start with two fields S Y and S /T' respectively on the surface of the sphere. As before we can 
take specific examples where these fields are either j± or We will keep the analysis generic here and will consider the specific examples later on. The 
spins associated with various fields are denoted by lower case symbols, i.e. u and v. The product filed now can be expanded in terms of the spin harmonics 
of spin u + v. as was done in Eq .d 1 6t . Similarly decomposing the other set of product field we obtain [ S »T" (il) s "iT'" (€l)]i m - We next construct the the 
power spectrum associated with the trispectrum from these harmonics: 



c ! rr '' r " r "'(r 1 ,r 2 ) = ^ £ [ s r s ,rV(n)[ s »rVr"']? m (r2) (36) 

m— — L 

This particular type of power spectra associated with trispectra has been studied extensively in the literature in the context of CMB studies; especially to 
separate out effects due to lensing of CMB from other secondaries. This is one of the two degenerate power spectra associated with trispectrum. After 
going through very similar algebra outlined in the previous section we can express the Cf T ,r r in terms of the relevant trispectra which it is probing. 
The expression in the absence of any mask takes the following form: 



3 „r", s ,„r"', s ,„r"' ej±,v,T,g (37) 

The Gaussian contribution from the disconnected part of the trispectrum G needs to be subtracted to construct the estimator. In the presence of a 
completely general mask w(Cl) the pseudo-Cf s or PCLs will have to be modified to take into account the effect of mask. This involves computing the 
spherical harmonics of the masked field [r(f2)r'(£2)u;(fj)] and cross-correlating it against the harmonics of [V" (Cl)T" (Cl)w(Cl)]. 



~rr'r"r'" 1 » 1 

Cf rs r (ri,r 2 ) = 21 + 1 E I sr s' r ' Hi"i[s" r "s'" r '" w]im\ Wl ^ 21 + 1 ^ w imW*i m (38) 

m— — I m — — l 

The resulting PCLs are linear combination of their all-sky counterpart. The mixing matrix which encode the information about the mode mixing will 
depend on the power spectrum of the mask as well as the spins of all four associated fields. The mixing matrix Mu< expressed in terms of the Wigner's 
3j symbols take the following form: 



W„; =^I>' + l)^ + l)( s ; s , o -(s' + s')){s"ls- V+*'")W; s,s',s",s"'e Q ,l,±2,3; 



(39) 

The spin indices s, s',. . . take values for k, ±2 for y+,1 for F and 3 for G. The pseudo-CV s expressed as a linear combination of all-sky power spectra 
can now be expressed using the following relationship: 
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~ rr / r"v'" i \ Tkfsa 1 's",s'" ^,rr',r"r"' t \ , ArY . 
V 

For near all sky surveys and with proper binning the mixing matrix M;;/ can be made invertible. This provides an unique way to estimate all-sky 
C,, ' and the associated information contents about the trispectra. For a given theoretical prediction the all-sky power spectra C\, ' can be 

analytically computed. Knowing the detailed m odel of an experimental mask allows u s to compute the observed Cf r r r accurately. The results 
presented here generalizes the ones obtained in dMunshi et aill2009l : iMunshi et alll201Cl) for the case of shear and Flexions. While the shear 7 like the 
polarization fields Q ±iU are spin -2 objects, flexions on the other hand can be associated with higher (or lower) spins. 

These results assumes generic field variables s T(r), S /F'(r) which can have arbitrary spin associated to them. We next specifiy certain specific cases 
where we identify three of the fields S T, s iT', s iiT" = * and s >i>V" = "f±. The other combinations can also be obtained in a similar manner. 



C % -2^ JhhiJhUiT 9hEli {I) ^ )[±2 +2 )' ( } 



4n ^ Jy ' I V ±2 =F2 

la. V ' V 

Other estimators for mixed trispectra involving different combinations of iJ-polarization and a tracer field (scalar) ^ can be derived in a similar manner 
and can provide independent information of corresponding trispectra. We have ignored the presence of B-type polarization in our analysis. Presence of 
non-zero B-mode can be dealt with very easily in our framework but the resulting expressions will be more complicated. 

3.2.2 The Three-to-One Power spectrum: C\ ' 

The number of power spectra that can be associated with a given multispectra depends on number of different way the order of the multispectra can be 
decomposed into a pair of integers. The bispectrum being of order three can be decomposed uniquely 3 = 2 + 1 and has only one associated power 
spectra. On the other hand the trispectrum whose order (= 4) can be decomposed in a two di fferent way; i.e., 4 — - 3+1 = 2 + 2. Hence trispectrum of a 
specific type correspond to a pair of two different power spectrum associated with it (see e.g. lMunshi et alj j2009) for more details). The results presented 
here are generalization for the case of non-zero spins. 

The other power spectrum associated with the trispectrum is constructed by cross-correlating product of three different fields [ s r s /r' s //r"] with the 
remaining field s mY"' . The cross correlation power spectrum in terms of the multipoles are given by the following expression: 

[ s r(r) s ,r'(r) s »r"(r)] im = J[ a T s ,V' a „T"\[ a+a , +s , l Y lm {Cl)*]dCl; [ a „,r'"(r)] im = [[ s »,r'%,,,Yi m (fl)*]d£l; (43) 
q rs ' r ' s " r "' s "' r "Vi,r 2 ) = _1_ [«r s) r' s „r"] im (ri)[ sW r'"]L(r 2 ). (44) 

m — — l 

By repeated use of the expressions Eq.([5j or equivalently Eq.lJSJl to simplify the harmonics of the product field in terms of the individual harmonics, we 
can express the all-sky result in the following form: 



_,rT'r",r'"/ \ r i ^^ ^ ^l ^ ^2/r^ f h h L \ f L h I 

H 7^2 ,^3 ,L 

s,s',s",s" € 0,1, ±2, 3. (45) 
In case of of partial sky coverage with a generic mask w(Cl) the relevant expression for the Pseudo-CV s will involve the harmonic transform of the mask. 



[ a T(r) s ,T'(r) s „T"(r)w(fl)]im = J [ S T S ,T' ^V'w^u+v+AiQ^dfl; 

[ s ,„T"'(r)w(Cl)]i m = J [ s -r"'H[ s ,»y" ; * m (fi)]dn; c<; r ,' r ' s " r "',"' r "' = M w cf°' T ' °" T " •°"' T "' (46) 

The mixing matrix has the following expression in terms of various spins involved and the power spectra of the mask introduced before. Note that the 
mixing matrix is of different form compared to what we obtained for two-to-one power spectra. This is related to how various fields with different spins 
were combined to construct these two estimators. In case of convergence trispectrum (spin-0) the two mixing matrices take the same form. 



^ V ' y '' = ^B^ + i)(^ + 1 )( (s + J + s «) o - (s + !'+,''))(*'" o i») |u ' ij2; (4?) 

The expressions derived here are derived for general mask but Gaussian noise of intrinsic ellipticity distribution of Galaxies. Any residual non- 
Gaussianity from systematics will have to be subtracted out. The estimators derived above indeed are unbiased but unoptimized. Optimization of these 
estimators will however have to be done in a model dependent way and from pre-constructed E and B maps from observed shear data. 



For a specific example we choose S T — S /T' — S "T" = 9 and S «/P 
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— E. In this case the three-to-one estimator takes the following form: 
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„***, 7± j T . h h L\(L fa I' 



^ = i£(*+w+i)(5 l ; l ; ilkft (48) 



47r " MO / V ±2 =F2 

The other possibilities include the pure trispectrum that corresponds only to _B-type polarization which can be used directly from a cubic shear map 
cross-correlated against a shear map: 



c i - 7 j JHl-ihJhHi 1 E h E l ( L ) I ±2 



^i fa L \ ( L fa I' 



±2 =p4 / V =F4 ±4 



4tt ^ v " ' V ±6 =F6 J V ±2 =F2 

The generalization to trispectra involving the flexions T and Q will follow a procedure similar to that of bispectrum. 

The bispectrum is defined through a triangular configuration in the multipole space. The trispectrum on the other hand is associated with a quadrangle 
in the multiple space. The extraction of information by the associated power spectra defined above corresponds to summing over all possible configurations 
keeping one of the sides of the triangle fixed. In case of power spectra associated with trispectra there are two different options: to keep the diagonal fixed 
and sum over all possible configurations C,' 2 ' 2 ' or to keep one of the sides fixed and sum over all possible configurations C;' 3 ' 1 '. 

It is possible to however introduce a window optimized to selectively search for information for a specific configuration either for bispectrum or for 
trispectrum. However such analysis which introduces mode-mode coupling can not be generalized to the arbitrary partial sky coverage as there is already 
a coupling of modes because of non-uniform coverage of the sky. 

Unlike the bispectra, the trispectra has non-vanishing contribution from the Gaussian component. The Gaussian contribution represents the un- 
connected component of the total Trispectra and needs to be subtracted out. These can be treated by using a Monte-Carlo approach which uses 
identical mask. In our analysis we have used the harmonic space approach which is complementary to real space analysis. It is possible to device 
statistics such as collapsed four-point functions, (~f 2 (Cl)^ 2 (Cl')) c or (~f 3 (Cl)'y(Cl)) c wh ich will generalize the bispectrum related statistics proposed by 
l lBerna rdeau. Me llier & Van Waer beke 2002; Bernarde au. van Waerbeke & M ellier 2003|) in real space, for the case of three-point correlation functions, 
to four-point correlation function in real space or equivalently to power spectrum related to trispectra. However higher order correlation function as well 
as their Fourier transforms will be more dominated by noise and higher number density of galaxies will be required to probe non-Gaussianity beyond the 
bispectrum. 

We also want to point out that at higher redshift, the number density of sources will be low. This will result in increase of shot noise. It is possible 
to include a redshift dependent weight function w(r). In the simplest case, this would reflect the number density of source galaxies. Such a weighting 
scheme can be incorporated in the modeling of our mask, which we have left completely arbitrary. If we use two different set of mask for two different 
redshift bands we can still use the expressions derived above for coupling matrices M U r. However we need to replace the power spectrum \wi | 2 of the 
mask with a cross-spectra of two different masks that are being used. 

It is possible to compute the covariance of our estimates of various Cjs, i.e. CiCy under certain simplifying assumptions which involves modeling of 
higher order correlations accurately. The shot noise contributions, involving various powers of source galaxy density, will dominate at higher I. A detailed 
analysis will be presented elsewhere. 



4 MODELLING GRAVITY INDUCED NON-GAUSSIANITY 

We start by modeling the bispectrum and trispectrum of the convergence field k(Ci) in terms of the underlying mass distribution 5. Next we will link the 
bispectrum and trispectrum associated with the shear "/(Cl) or flexions J-(Cl) or Q(Cl) fields in terms of the corresponding quantities for the convergence 
7. These will be used for the computation of the power spectra associated with these multispectra of various order. 



4.1 Linking Weak Lensing Observables and Mass Correlation Hierarchy 

The two-point and higher order correlation functions as well as their Fourier counterparts, the power spectrum, the bispec trum and other higher order 
multispectra beyond the lowest orders are well understood for weak lensing convergence field jMunshi & Jair]|200dl200lh . The modeling depends on 
accurate modeling of underlying mass distribution. The projected convergence power spectrum Cf can be expressed in terms of the 3D mass power 
spectrum P s : 

Here <1a (r) is the angular diameter distance at a comoving distance r, corresponding scale factor is a and r s is the distance to the background source 
distribution (we assume that all sources are situated at a redshift of z s = 1). The power spectrum of the convergence field k is therefore a weighted 
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projection of the 3D matter power spectrum along the line of sight. The weight w(r) depend on underlying cosmology through the Hubble constant Ho, 
!)m as well as cLa- It also depends on the source redshift. 



(Ki imi Ki 2 m 2 Ki 3 m 3 )c - ( mi ma ma J B h i 2 i 3 ; B hh i 3 — ( mi ma rn 3 J (^ii™!^™?^^)? (51) 

m\m 2 m^ 

The bispectrum of the convergence can now be expressed as the bispectrum of the underlying mass distribution. The all-sky expression for bispectrum 
can be written as: 



RKKK _ r f,y(r) R if h h h \ / (2? 1 + l)(2/ 2 + l)(2/3 + l) f h k h\ 

A similar expression holds for trispectrum Xj Yi (L) in terms of matter trispectrum. The above expressions are derived within the context of Limber's 
approximation as well as the Born's approximations which are widely used for such calculations. The flat sky bispectrum bi t i 2 i 3 is related to its all-sky 
counterpart by the relation Bi t i 2 i 3 = Iixl^la^ixinh' Using expression derived previously Eg. d 15b or Eq.J 1 3t we can relate the bispectrum for convergence 
and the bispectrum associated with the shear or flexion field are related B^^f 3 = F^F^F^Bf^ la with F ; B = J ^"^^ y^-. For large values of 
I Fi ~ 1 and different power spectra that we defined associated with bispectra differ only through the geometrical form factor as encoded in the Wigner 
3j symbols. 



4.2 Matter Correlation Hierarchy 

It is clear that we need accurate analytical modeling of dark matter clustering for prediction of weak lensing statistics. We lack detailed such theoretical 
understanding of gravitational clustering. On larger scales, where the density field is only weakly nonlinear, perturbative treatments are kn own to be valid. 
For a phenomenological statistical description of dark matter clustering in collapsed objects on nonlinear scales, typically the halo model dCoorav & Sethi 
2002) is used. We will be using the halo model in our study. However, an alternative approach on small scales is to employ various ansatze which trace their 
origin to field theoretic techniques used to probe gravitational clustering. The hierarchical ansatz has also been used for many weak lensing related work, 
where the higher-order correlation functions are constructed from the two-point correlation functions. Assuming a tree model for the matter correlation 
hierarchy (typically used in the highly non-linear regime ) one can write the mo st general case, the N point correlation function, (f N {r\, . . . , r„) as a 
product of two-point correlation functions fl(| r i — r j|) dBernardeau et all l2002V Equivalently in the Fourier domain the multispectra can be written as 
products of the matter power spectrum Ps(ki). The temporal dependence is implicit here. 

(iV-l) 

^(r 1 ,...,r ) = ({(r 1 )..i(r n )) c = £ Q N , a ]T ]J fc(|r 4 - rj -|). (53) 

ck,N— trees labcllings cdgcs(ij) 

It is however very interesting to note that a similar hierarchy develops in the quasi-linear regime at tree-level in the limiting case of vanishing variance. 
However the hierarchical amplitudes become shape-dependent in such a case. These kernels are also used to relate the halo-halo correlation hierarchy 
with underlying mass correlation hierarchy. Nevertheless there are indications from numerical sim ulations that these amplitudes become co nfiguration- 
independent again as has been shown by high resolution studies for the lowest order case Q3 = Q dScoccimarro et "a ll998l ; Bernardeau et aill2002h . bee 
Waer bekeetaj d200lh for related discussion about use of perturbation theory results in intermediate scales. In the Fourier space however such an ansatz 
means that the entire hierarchy of the multi-spectra can be written in terms of sums of products of power spectra with different amplitudes Qn,c etc. The 
power spectra is defined through (5(ki)5(k2)} e = (2-7r) 3 6>3.D(ki 2 )P(fci). Similarly the bispectrum and trispectrum are defined through the following 
expressions (<5(ki)<5(k 2 )<5(k 3 )) c = (2 7 r) 3 5 3 i3(ki 2 3)B' 5 (ki, k 2 , k 3 ) and (<5(ki) ■ ■ ■ <5(k 4 )) c = (27r) 3 <33i ? (ki 2 34)T i (ki,k 2 ,k3,k4). The subscript c 
here represents the connected part of the spectra and k^ ...^ = k^ + • • • + kj ra . The Dirac delta functions S^d ensure the conservation of momentum at 
each vertex representing the multispectrum. 



B 4 (k 1 ,k 2 ,k 3 ) Ekl= o = Q 3 [P S (ki)P S (k2) + P S (k 1 )P s (k 3 ) + P s (h 2 )P s (k 3 )] (54) 
T\^ 1 MMM)Y.^=o^Ra[P\k l )P\k 2 )P\k 3 ) + c y c. P erm] + R b [P\k 1 )P\\^ 1 2\)P\ (55) 

Different hierarchical models differ in the way numerical values are allotted to various amplitudes. iBernardeau & Schaeffed dl992l) considered "snake", 
"hybrid" and "star" diagrams with differing amplitudes at various order. A new "star" appears at each order, higher-order "snakes" o r "hybrid" diagrams 



are built from lower-order "star" diagrams. In models where we only have only star diagrams dValageas. Barber. & Munsh i 2004) th e expressions for 
the tri spectrum takes the following form: T s (ki, k 2 , k3, krt)^ k .= = Qa\P s (^ ) P s (k?) P s ( kg) + cvc.verm.]. Following lValageas. Barber. & M unshi 
d2004l) we will call these models "stellar models". Indeed it is also possible to use perturbative calculations which are however valid only at large scales. 
While we still do not have an exact description of the non-linear clustering of a self-gravitating medium in a cosmological scenario, theses approaches 
do capture some of the salient features of gravitati onal clustering in the highly non-li near regime and have been tested extensively against numerical 
simulation in 2D statist ics of convergence of shea r dValageas. Barber. & Munshi 2004). These models were also used in modeling of the covariance of 
lower-order cumulants dMunshi & V alageas 2005]). 
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4.3 Halo Model 

Halo model on the other hand relies on modeling the clustering of halos and predictions from perturbative calculation to model the non-linear correlation 
functions. It is known that t he halo over density at a given position x, 5 (x, M ; z) can be related to the underlying density contrast <5(x, z) by a Taylor 
expansion as was shown by dlVlo. Jing & White! 19971) . 

S h (x, M; z) = h (M; z)<S(x, z) + ifo 2 (M, z)S 2 (x, z) + . . . (56) 

The expansion coefficients are functions of the threshold v c — S c /a(M, z). Here S c is the threshold for a spherical over-density to collapse and <r(M, z) 
is r.m.s fluctuation within a tophat filter. The halo model incorporates the perturbative aspects of gravitational dynamics by using it to model halo-halo 
correlation hierarchy. The nonlinear features take direct co ntribution fro m the halo profile. Other ingredients being number density of halo. The total 
power spectrum P t (k) at non-linear scale can be written as ( Seliak 2000) 

P PP = I°{k,k); P HH (k) = [ll(k)] 2 P(k); P* = P HH (k) + P pp {k). (57) 

In our calculation minimum halo mas that we consider is 10 3 M and the maximum is 110 16 M. It is known that the massive halos do not consider 
significantly due to their low abundance. The bispectrum involves terms from one, two or three halo contributions and the total can be written as: 

B* (fci , k 2 , k 3 ) = B ppp (ki , k 2 , k 3 ) + B ppp (fci , k 2 , k 3 ) + B ppp (ki , k 2 , k 3 ) ; (58) 
B PPP = Il{ki,k2,k 3 ); B PHH (k 1 ,k 2 ,k 3 ) = I 2 L (k 1 ,k 2 )I°(k 3 )P{k 3 ) + cyc.perm.; (59) 
B HHH (k 1 , fca, k 3 ) = [2J(fei, k 2 , k 3 )ll{k 3 ) + / 1 2 (fc3)]/ 1 1 (fc 1 )/ 1 1 (fc 2 )P(fc 1 )P(A :2 ) + cyc.perm. (60) 

The kernel J(ki, k 2 , k 3 ) is derived using second order perturbation theory [|ry| dl984h ; iBouchet et all dl992h which we discussed above in the context of 
hierar chical ansatz. The integrals 1^ can be expressed in terms of the Fourier transform of halo profile (assumed to be an NFW dNavarro. Frenk & White! 
1996) profile): 



/,/lA-t. <'••:-• A-,,:-.) = / ,IM ( ^ dn{ ^ z) b^M)y(k 1 ,M)...y(k^M); y(k, M) = JL £ dr4nr 2 p(r, M) 



\(kr) 



(61) 



At 

The mass function is assumed to be given by Press-Schechter mass function jPen et alf2 003). Our calculations are based on using the expressionsl60lin 
equation Eq.J50t and Eq.J52t, The convergence power spectra and bispectra are next related to the observables involving shear using Eq.j27}.Eq<22) and 
Eq.Olll. We will only be considering only lowest order non-Gaussianity statistics, i.e. bispectrum here but results can readily be generalized to higher 
orders. 



5 CONCLUSIONS 

Future weak lensing surveys will play a crucial role in cosmological research, particularly in further reducing the uncertainty in fundam ental properties of 
the standard cosmological model, including those that describe the evolution of equation of state of dark energy lRefregieretaI] ( l2010h . It is known that, 
by exploiting both the angular diameter distance and the growth of structure, wea k lensing surveys are able to const r ain cosmological parameters, as we ll 



as testing the accuracy that general relativity provides as a description of gravity Heavens, Kitching & Verdel d2007t) : lAmendola. Kunz & Sapond J20081) : 
iBenvon. Bacn & Kovamal d2009h : ISchrabback et all d2009h : lKiibinger e t alj d2009h . 

Constraints from weak lensing surveys are complementary to those obtained from cosmic microwave background studies and from galaxy surveys 
as they probe structure formation in the dark sector at a relatively low redshift range. Initial studies in weak lensing were restricted to studying two-point 
functions i n projectio n for the entire source distribution. It was, however, found that binning sources in a few photometric redshift bins can imp rove the 



tunctions in projection tor the entire source distribution, it was, however, touna that binning sources in a tew photometric redshitt bins can improve the 
constr a ints iHul (1 1 999 ). Mor e recently a full 3D fo rmalism has been developed which uses photometric redshifts of all sources without any binnin tjHeavensI 



d2003l) :l Castro et aif l2005l) :l Heav ens et all d2006h . These studies h ave demonstrated that 3D lensing can provide more powerful and tighter constraints 



on the dark energy equation of state parameter, on neutrino masses Ide Bernardis et al.1 J2009T) . as well as testing braneworld and other alternative gravity 
models. Mos t of these 3D works have primarily focu s ed on power spectrum analysis, but in future accurate higher-order statistic measurement should be 
possible (e.g. lTakada & Jainl j2004):ISemboloni et all fe009h ). 



In a recent work ( Mun shi. Heavens & Colesll2010l)~ extended previous studies using convergence maps to probe gravity-induced non-Gaussiniaty. 
These studies go beyond previous analysis of power-spectrum estimation from convergence data. Recovery of cosmological information however is even 
more straightforward from shear maps. In this paper we have generalized the results previously obtained for convergence to shear. We have gone beyond 
conventional power spectrum analysis by defining power spectra associated with higher-order multispectra in such a way as to compress information 
available in the multispectra to a single derived power spectrum. Generic results which do not depend on detailed modeling of multispetcra have been 
obtained and later specialized with concrete models using specific forms for the gravity-induced clustering as encoded in the higher-order multispectra. 
The generic results we have presented do not depend on the specific nature of non-Gaussianity and are able to handle primordial as well as secondary 
non-Gaussianity. 

Though the higher-order multispectra contain a wealth of information in through their shape dependence, they are partly degenerate. It would be 
desirable to exploit the entire information content encoded in higher order multispectra for constraining structure formation scenarios. However, the direct 
determination of the multispectra and their complete shape dependence from noisy data can be a very difficult t ask. In this paper we have emp loyed a set 
of statistics called "cumulant correlators" which were first used in real space in the context of galaxy surveys dSzapudi & Szalavll 1993lfT997h and later 
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extended to CMB studies dMunshi et all20ld.l2009h . We have presented a general formalism for the study of the power spectra or the Fourier transforms 
of these correlators. We present a 3D analysis which takes into account the radial as well as on the surface of the sky decomposition. Previous studies in 
weak lensing mainly concentrated on convergence maps containing spin-0 objects that can be treated relatively easily. However, from the point of view 
of observational data reduction, it is much more natural to focus on galaxy shear which is a natural byproduct of weak lensing surveys. In this paper we 
have extended the concept of cumulant correlators to higher spin objects to shear as well as flexions. We consider the correlations of cumulants which 
are constructed from such spinorial objects. The analytical results we obtain are valid for an arbitrary mask and our formalism allows us to correct for 
the effect of mask at arbitrary order for general spinorial field. The analytical results are sufficiently generic to include a tracer field for the large scale 
structure and probe mixed bispectrum involving underlying mass distribution as well as the large scale tracer field. 

We have restricted this study to the third and fourth order, as data from observations gets increasingly noise dominated as the order increases. 
However, our analysis can be generalized to higher order and some of our results are indeed valid at arbitrary order. At third order, we define a 
power spectrum which compresses information associated with a bispectrum to a power spectrum. This power spectrum Cf F ' r (r2,r\) is the cross- 
power spectrum associated with squared convergence maps r(ri, Q)r'(ri, Cl) constructed at a specific radial distance r\ against T"(r2,Cl) at V2- In 
a similar manner we also associate power spectra C l ' and C l ' with associated trispectra (L; ri). There are two different power 
spectra at the level of trispectra which are related to the respective real-space correlation functions (r(n, Q)V"(ri, Cl)F"(r2, Q')V"(r2, Cl)) and 
(r(ri, f2)r'(ri, Cl)T"(ri, Cl)F'"(r2, Cl')). We expressed these real-space correlators in terms of their Fourier space analogue which take the form of 
Cf r ' r "' r "'(r 2! r 1 )andCf r '' r " r "'(r 2) r 1 ). 

Th e statistical descrip tors we have presented here will provide particularly useful tools for the study of non-Gaussianity in alternative theories of 
gravity (Bernardeau 2004); which is one of the important science drivers for the future generations of weak lensing surveys. We plan to present detailed 
results elsewhere in future. The various analytical ap proximations, such as the Born approximations linking the observed shear or convergence bispectrum 
(or trispectrum) have been studied with some rigor dHamana et a l. 2002; Shapiro & Coorav 2006). The effect of source clustering, as well as lens and 
source overlap - which also introduce corrections - have also been studied with some detail. We plan to extend these studies for the power spectra presented 
here elsewhere. 

The estimators we have worked with can be generalized further to take into account optimum weighting. This will involve inverse covariance 
weighting the observed field harmonics. This will then be useful in constructing a matched filtering version of the estimator used here to fine tune 
the study of primordial as well as gravity induced non-Gaussianity. While previous ray tracing simulatio n of weak lensing used a flat-sky approach 
jjain. Seliak & Whitdb oOO). simulations are now available for the entire observed sky dTevssier et al.ll2009T) . This will provide an opportunity to test the 
analytical results presented here. A detailed discussion will be presented elsewhere (Munshi et al. 2010, in preparation). 
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